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Is There a Missing Galaxy Problem at High Redshift? 

Kentaro Nagamine^, Renyue Cen^, Lars Hernquist^, Jeremiah P. Ostriker^'^, 

& Volker Springel^ 

§ ■ ABSTRACT 
(N 

We study the evolution of the global stellar mass density in a A cold dark 
matter (ACDM) universe using two different types of hydrodynamical simulations 
(Eulerian TVD and SPH) and the analytical model of Hernquist & Springel 
(2003). We find that the theoretical calculations all predict both a higher stellar 
^ mass density at 2 ~ 3 than indicated by current observations, and that the 

>- ■ peak of the cosmic star formation rate history should lie at 2 > 5. Such a star 

^ ■ formation history implies that as much as (70%, 30%) of the total stellar mass 

CN ■ density today must already have formed hj z = (1,3). Our results suggest that 

current observations at z ~ 3 are missing as much as 50% of the total stellar 
. mass density in the Universe, perhaps owing to an inadequate allowance for dust 

. obscuration in star-forming galaxies, limited sample sizes, or cosmic variance. 

p We also compare our results with some of the updated semi-analytic models of 

Q . galaxy formation. 

^ ■ Subject headings: cosmology: theory — stars: formation — galaxies: formation 



galaxies: evolution — methods: numerical 



1. Introduction 



Is the evidence for high redshift galaxy formation consistent with the concordance 
ACDM model? Recent observational results include the discovery of Extremely Red Ob- 
jects at z>l (e.g. Elston et al. 1988; McCarthy et al. 1992; Hu & Ridgway 1994; Cimatti 
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et al. 2003; Smail et al. 2002), Sub-millimeter galaxies ect z > 2 (e.g. Small et al. 1997; 

Chapman et al. 2003), Lyman-break galaxies (LBGs) at 2; ~ 3 (e.g. Stcidcl ct al. 1999), and 
galaxies at ^ > 4 either by their Lyman-a emission (e.g. Hu et al. 1999; Rhoads & Malhotra 
2001; Taniguchi et al. 2003; Kodaira ct al. 2003; Ouchi et al. 2003a) or by their optical to 
near infrared (IR) colors (e.g. Iwata et al. 2003; Ouchi et al. 2003b; Dickinson et al. 2003b). 

Multiband photometry including the ncar-IR band makes it possible to estimate the 
stellar mass of these high redshift galaxies by fitting the observed photometric results with 
artificial galaxy spectra generated by a population synthesis model. Using this technique, 
several groups have now estimated the stellar mass density in the Universe in the redshift 
range of < ;2 < 3 (e.g. Brinchmann & Ellis 2000; Cole et al. 2001; Cohen 2002; Dickinson 
et al. 2003a; Fontana et al. 2003; Glazebrook et al. 2004). Rudnick et al. (2003) also estimated 
the stellar mass density at 2; = — 3 by combining the estimates of the rest-frame optical 
luminosity density and the mean cosmic mass-to-light ratio. These observational estimates 
constrain the evolution of the stellar mass density ft^ as a function of redshift or cosmic 
time. 

The observations indicate the presence of a significant stellar population at high redshift, 
and, by comparing with semi-analytic models of galaxy formation, some authors claim that 
ACDM models seriously underpredict galaxy formation at z ~ 3. For example, Fontana et al. 
(2003) compare their estimate of Q^, from the Hubble Deep Field (HDF) South with the semi- 
analytic model of Mcnci ct al. (2002), and argue that the high-mass tail of the galaxy stellar 
mass function is not adequately described by CDM models. Comparing to the same semi- 
analytic model. Poll et al. (2003) argue that hierarchical models lack sufficient star formation 
aiz — 2 — 4, resulting in a failure to reproduce the pronounced brightening of the luminosity 
function at these redshifts. Also, Dickinson et al. (2003a) find that their data from the HDF 
North suggest a steeper increase in fi^ than some semi-analytic models (Kauffmann et al. 
1999; Somerville et al. 2001; Cole et al. 2000), and some of the semi-analytic models predict 
higher stellar mass density compared to their data points. Given these contradicting claims 
and the large parameter space available to the semi-analytic models as well as the limitations 
of the current observational samples, it is not clear at present if any of these discrepancies 
actually pose a serious problem to hierarchical evolution. 

In this paper, we compare the observational data with the results from state-of-the- 
art cosmological simulations of the standard concordance ACDM model and the theoretical 
model of Hernquist & Springel (2003) (hereafter H&S model), to show that, contrary to some 
of the claims, theory predicts a higher fl-^ at z = 3 than indicated by current observations 
and that the cosmic star formation rate (SFR) density peaks at 2; > 5, earlier than sug- 
gested by most semi-analytic models. Also, we will explicitly compare our results with those 
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from the updated semi-analytic models by Somerville et al. (2001), Granato et al. (2000) 
(GALFORM), and Menci et al. (2002). 

2. Simulations 

We will show results from two different types of cosmological hydrodynamic simulations. 
Both approaches include "standard" physics such as radiative coohng/heating, star forma- 
tion, and supernova (SN) feedback, although the details of the models and the parameter 
choices differ somewhat. 

One set of simulations was performed using an Eulerian approach, which relies on a 
particle-mesh method for the gravity and the Total Variation Diminishing (TVD) method 
(Ryu et al. 1993) with a fixed mesh for the hydrodynamics. The treatment of the radiative 
cooling and heating is described in Cen (1992) in detail. The structure of the code is similar 

to that of Cen & Ostrikcr (1992, 1993), but the code has significantly improved over the 
years with additional input physics. It has been used for a variety of studies, including the 
evolution of the intergalactic medium (Cen et al. 1994; Cen & Ostrikcr 1999a,b), damped 
Lyman-a absorbers (Cen et al. 2003), and galaxy formation (e.g. Cen & Ostriker 2000; 
Nagamine, Fukugita, Cen, & Ostriker 2001a,b; Nagamine 2002). 

Our other simulations were done using the Lagrangian Smoothed Particle Hydrodynam- 
ics (SPH) technique. We use an updated version of GADGET (Springel et al. 2001), which 
uses an 'entropy conserving' formulation (Springel & Hernquist 2002) to mitigate problems 
with energy /entropy conservation (e.g. Hernquist 1993) and overcooling. This code also uses 
a multiphase model of the interstellar medium to describe self-regulated star formation and 
a phenomenological model for galactic winds (Springel & Hernquist 2003a). This approach 
has been used to study the evolution of the cosmic SFR (Springel & Hernquist 2003b), 
damped Lyman-a absorbers (Nagamine, Springel, & Hernquist 2004a,b), and galaxies at 
high redshifts (Nagamine, Springel, Hernquist, & Machacek 2004). 

The cosmological parameters adopted in the simulations are intended to be consistent 
with recent observational determinations (e.g. Spergel et al. 2003), as summarized in Table 1. 

3. Stellar Mass Density 

Figure 1 shows the evolution of the global stellar mass density Q^, as a function of 
redshift. Figure la compares the simulation results with the observations, and Figure lb 
compares the semi-analytic models with the observations. Our results are scaled to unity at 
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z — (see discussion below) . Observational data points are normalized to the local estimate 
by Cole et al. (2001, filled circle at z=0) following Rudnick et al. (2003). The result of the 
H&S model shown in both panels is simply the integral of their approximate formula for the 
SFR which will be given in Equation (1). 

For recent epochs {z <1), the computed results are within the range of the observational 
data. However, the important result here is that at 1 < z < 3, all the observational estimates 
(both corrected and uncorrected estimates for the incompleteness of the survey) are smaller 
than the simulation results by more than a factor of two. The early development of the stellar 
mass density is much faster in the simulations than is suggested by current observations. As 
we discuss in Section 4, the larger fi^ at 2; ~ 3 originates from a higher SFR at z > 3 which 
peaks at 2; > 5 in both the TVD and SPH simulations. Such a star formation history implies 
that about 70% (50, 30, 15%) of the total stellar mass density today must have been in place 
hy z — 1 {z — 2,3,5). The contribution to Qi,{z — 0) from the star formation at 2; > 6 is 
about 10%. 

We note that the data points of Rudnick et al. (2003, magenta filled squares) should be 
considered as a lower limit, as they did not attempt to make any extrapolations to correct for 
the incompleteness of their data. On the other hand, Dickinson et al. (2003a, green square 

boxes) correct for the incompleteness of their data by integrating the Schechter fit with a 
faint-end slope of a = —1.4 down to a fainter magnitude. For the Fontana et al. (2003, 
blue crosses) data points, we have plotted the values for the SMC extinction case (which 
yields slightly larger values compared to the Calzctti ct al. (2000) extinction), and the 
upper error bars are extended up to their 'maximum mass' estimates. For the Glazebrook 
et al. (2004, black filled triangles) data points, we plotted the values for the mass limit of 
log(M/M0) = 10.2 which yields largest estimates from their data. We also note that dust 
extinction is taken into account in the analyses by Brinchmann & Ellis (2000); Dickinson 
et al. (2003a); Fontana et al. (2003); Glazebrook et al. (2004) by allowing the extinction 
parameter to vary when fitting the broadband photometric measurements of each galaxy 
by the spectral energy distribution generated by a population synthesis model. Rudnick 
et al. (2003) adopt the Calzetti extinction with E{B — V) = 0.30 when deriving the mean 
mass-to-light ratio from the rest-frame U — V color. 

The scaling of the simulation results to the local estimate at ^ = leaves us with some 
concerns, because the relative value of docs indeed depend on this scaling. However, some 
kind of normalization is necessary to compare the results of different numerical models on the 
same basis, because every numerical model has its input parameters, and what we would like 
to focus in this paper is the relative speed of the development of Q^, from high redshift to the 
present time. The two TVD simulations in fact yield different values of fi* owing to cosmic 
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variance and differences in the set of cosmological parameters adopted in the simulations: 

= 0.0077 and 0.0052 for N864L22 and N768L25, respectively. The normalization of 
at z = in Figure 1 can be partially explained by the scaling of with the baryon mass 
density in the different simulations. Hernquist & Springel (2003) have shown that scales 
as ~ Q^'^ from theoretical arguments, and Gardner et al. (2003) found that the amount 
of cold gas and stars in their SPH simulations follow oc 0^ °^^ °. 

Using the cosmological parameters of the two TVD runs, the scaling fi^ oc Q^'^ gives 
the expected ratio r2^(N864L22)/f2^(N768L25) = 1.7. The actual corresponding ratio from 
the two simulations is (0.0077/0.0052) = 1.5, which is a reasonable agreement, given other 
uncertainties, such as cosmic variance. It is expected that the H&S model (blue long-dashed 
line in Figure 1) should yield the highest D.^, at z — 3, because it is intended to remove the 
effects of hmited resolution and cosmic variance. Since the results of N864L22 and N768L25 
are not corrected for the limited box-size and resolution, it is natural that they lack the 
earliest star formation at z > 10, resulting in a lower fl^ at z = 3. This is also true of 
the SPH G6 run which by itself cannot resolve the entire starforming population at z > 3. 
Also, a comparison of the N864L22 and N768L25 results gives an idea of the level of cosmic 
variance for a volume of ^ {25h^^ Mpc)^, and the deviation of the SPH G6 run from the 
H&S model is also a consequence of cosmic variance. 

Figure lb compares the results of H&S theoretical model, semi-analytic models of 
Somcrvillc ct al. (2001), Granato et al. (2000), Menci et al. (2002), and observations. The 
model shown for Somerville et al. (2001) is the updated 'accelerated quiescent' model which 
was also used in Somerville et al. (2004), and follows the merger tree down to halos with 
circular velocity 30 km s~^. The result shown for Menci et al. (2002) is the original model 
used by Poll et al. (2003) and Fontana et al. (2003) for the comparison to their observational 
data. Recently Menci et al. (2003) updated their model to include a starburst mode of star 
formation, which gives a result close to that of GALFORM. The results of the semi-analytic 
models are somewhat higher than the observational estimates, but within the upper limit of 
Fontana et al. (2003) data point at z = 3. It is also clear that the predicted by the H&S 
model is higher than that of the two semi-analytic models by about a factor of two. The 
discrepancy between the models and observations increases at ^ ~ 1.5, and is significant at 
more than 1-a. 



4. Cosmic Steir Formation Rate at < 2; < 6 

In Figure 2, we show the cosmic SFR density as a function of redshift. Panel (a) 
shows the simulation results, and panel (b) shows the semi-analytic model results. The 
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simulation results shown in this figure are extracted directly from the runs, without any 
further adjustments. The H&S model that is shown in both panels has an approximate form 
as ^ 

p* = p^(0)— — ^ . (1) 

where x{z) = {H{z)/ HqY/^. For a ACDM universe with the star formation and feedback 
algorithm described by Springel & Hernquist (2003a), the parameters defining the SFR 
density take the values a = 0.012, (5 = 0.041, and p*(0) = 0.013 M© yr"^ Mpc^^ See 
Section 5 for further discussion on this formula. 

The line types and data points are described in the caption. The observationally es- 
timated ultra-violet (UV) luminosity densities puv have been converted into the SFR by 
puv[erg s-^ Hz-^ Mpc"^] = 8.0 x IQ'^'^SFR[Mq yr"^ Mpc~^] (Madau, Pozzetti, & Dickinson 
1998). The SFR data points arc corrected for dust extinction according to a prescription 
similar to that of Steidel et al. (1999): we assume the highly uncertain extinction correction 
factors to be 1.3 {z < 2) and 4.0 (^ > 2), while Steidel et al. (1999) used higher values 2.7 
{z < 2) and 4.7 {z > 2). 

It is important to note that each data point has been derived under different assump- 
tions, because the faint-end slope of the luminosity function of galaxies at high rcdshift is not 
well constrained and adopting a steeper slope and integrating down to fainter magnitudes 
would certainly yield a larger UV luminosity density. Here we describe some of the high 
redshift data points. Steidel et al. (1999, open stars a.t z = 3, 4) derived their UV luminosity 
density by integrating the luminosity function with a faint-end slope of a — —1.6 down 
to O.IL*. For their data points, we read off the SFR from Fig.9, corrected to our fiat-A 
cosmology, and apphed our extinction correction. Giavahsco et al. (2004, open triangles at 
^ = 3 — 6) integrated the Schechter fit with a faint-end slope of o; = —1.6 down to O.2L3, 
where L3 is the characteristic UV luminosity of LBGs at z ~ 3. Iwata et al. (2003, open 
pentagon at ^ = 5) integrated their Schechter fit with a = —1.5 in the magnitude range 
of 22.5 < Muv — 51og/i < —20.0. This integration range was chosen to match with the 
limiting magnitude of the z ^ 4 sample by Steidel et al. (1999) in terms of absolute mag- 
nitude, therefore it is expected from the shape of the luminosity function that the resultant 
UV luminosity density at 2; ~ 5 is 0.56 — 0.69 times of that at 2; ~ 3 in the same absolute 
magnitude range, depending on the choice of cosmology and the integration range. For the 
Ouchi et al. (2003b, inverted open triangles) data points, we show their total UV luminosity 
density case where the luminosity function with a faint-end slope of a = —2.2 is integrated 
down to O.IL*. 

An important theoretical result here is that all the simulation results and the H&S 
model peak at 2; > 5, and not at lower redshifts, as is often found in semi-analytic models of 
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galaxy formation (e.g. Baugh et al. 1998; Kauffmann et al. 1999; Cole et al. 2000; Somerville 
et al. 2001, see panel (b)). Observationally, it is an unsettled problem whether the SFR levels 
off at high-redshift (e.g. Steidcl ct al. 1999) or still increases beyond z > 3 (Lanzetta et al. 
2002). In particular, Lanzetta ct al. (2002) stress the importance of the cosmological surface 
dimming effect and argue that previous measurements have missed a significant fraction of 
the ultraviolet luminosity density of the universe ai z >2. The scatter in the data points at 
2; > 3 is still large, and it is not possible at this time to determine the trend in SFR at these 
redshifts. However, our point is a theoretical one, and does not depend on the distribution 
of data points as we detail in Section 5. 

In passing, we note that the absolute values of at 2; = computed by integrating 
the SFR curves shown in Fig.2b are 0.0037, 0.0031, 0.0041, and 0.0090 for the models by 
Hernquist & Springel (2003), Somerville et al. (2001), Granato et al. (2000), and Menci et al. 
(2002), respectively. 

The result of TVD N864L22 at z < 1 is shghtly higher than most of the observational 
data points and has two large bumps at 2; ~ 1.0 and 2.0. These two peaks presumably 
owe to major merger events taking place in the simulation, and are a consequence of cosmic 
variance in the relatively small box utilized. If we had utilized a larger box, then we would 
expect the curve to be smooth. Clearly the box-size of Lbox = 20 — 30h~^ Mpc is not large 
enough to accurately model the volume averaged quantities in the universe at z < 1. 

5. Discussion & Conclusions 

We have shown that two independent different types of numerical hydrodynamic simu- 
lations both predict that the cosmic star formation rate density should peak at z > 5, and 
that this relatively early peak in the SFR leads to a more rapid development of the stellar 
mass density than current observational estimates. When all the results are scaled to the 
local value at z = 0, the stellar mass density at 2; = 3 in the simulations is larger than 
observed values by more than a factor of two. We also showed in Figure lb that the semi- 
analytic models of Somerville et al. (2001), Granato et al. (2000), and Menci et al. (2002) 
predict larger than current observational estimates, but within the upper error bar of the 
Fontana et al. (2003) data point at 2; = 3. In particular, the H&S model predicts higher 
than those semi-analytic models by about a factor of two at 2 = 3. This comes from the 
fact that both the simulations and H&S model has a peak of SFR at 2 = 5 — 7, whereas the 
semi-analytic models have a peak at z — 2 — A. 

The high predicted by our simulations and the model of Hernquist & Springel (2003) 
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suggest that current observations are missing nearly half of the total stellar mass density in 
the universe at high redshifts. This missing stellar mass could be hidden in a population 
of red galaxies that have not been detected previously in the optical ground-based data. In 
fact, Franx et al. (2003) and Daddi et al. (2004) find such a population of red galaxies at 
z > 2 in HDF-South that has a volume density half that of LBGs at 2; = 3 and a stellar mass 
density comparable to that of LBGs. However, the data points by Rudnick et al. (2003) 
include this red galaxy population. Therefore the cosmic variance might be a stronger cause 
for the underestimate of Jl^ by the current observations. This is hinted by the absence of this 
red population in HDF-North. If these red populations are strongly clustered as suggested 
by Daddi et al. (2003), observations with small fields-of-view could easily miss them. The 
stellar mass density in the simulations and observations agree reasonably well at 2; < 1 where 
the rate of increase is much slower than at higher redshift. In addition, the dust extinction 
correction may seriously underestimate the fraction of early star formation that is heavily 
obscured. If we only had optical UV observations of nearby starburst galaxies - rather 
than the full spectrum extending to the submillimeter - we would greatly underestimate 
star formation in these systems (e.g. Hughes et al. 1998; Meurer et al. 1999; Barger et al. 
2000; Calzetti et al. 2000; Takagi et al. 2003). Future studies of starforming galaxies in the 
far-infrared wavelengths (e.g. Kennicutt et al. 2003) using the Spitzer Space Telescope, in the 
submillimeter (e.g. Tecza et al. 2004), in the millimeter (e.g. Bekki & Shioya 2000) using the 
Atacama Large Millimeter Array (ALMA), and in X-rays (e.g. Norman et al. 2004) using the 
Chandra X-ray Observatory would enable us to constrain the nature of starburst galaxies 
better, and estimate the intrinsic SFR more accurately. 

It is reassuring that the two different sets of simulations, which use very different hy- 
drodynamic methods (i.e. Eulerian TVD and SPH), both give a star formation history that 
peaks at 2; > 5. As argued by Hernquist & Springel (2003) this is to be expected, because 
the evolution of the cosmic SFR is driven mainly by a competition between gravity and the 
expansion of the Universe, with a weaker dependence on the details of star formation and 
feedback. The form of equation (1) given in Section 4 can be understood as follows. At 
high redshifts, when the cooling time is short, star formation is limited primarily by the 
gravitational growth of halos, which is independent of the dissipative gas dynamics. Thus, 
the parameter /3 and the factor x^/^ in the exponential of equation (1) are determined by the 
form of the matter power spectrum, and the description of star formation and feedback enter 
into p-k{z) only logarithmically. At low redshifts, the supply of star forming gas is limited 
by the expansion rate of the Universe, fixing the dependence of Pi,{z) on x as 2; ^ 0. The 
explicit influence of the prescription for star formation and feedback is again subdominant 
and mainly affects the values of the normalization parameters a and p^(0). 

For these reasons, equation (1) can be generalized straightforwardly to other cosmologies 



-9- 



and to include other physics (e.g. Yoshida et al. 2003). Moreover, as shown in e.g. figure 
6 of Hernquist & Springel (2003), the fact that the SFR density is regulated mainly by the 
competition between gravity and the expansion of the Universe means that the peak in Pi,{z) 
should lie at 2; ^ 5, unless an implausible value is adopted for the parameters governing 
star formation. Hence, it is not surprising that the various sets of simulations should be 
consistent, with residual differences owing to details in the cosmology and, most important, 
cosmic variance and resolution hmitations. We plan to investigate these issues in the future 
using the algorithms described here, as well as adaptive mesh refinement codes. 

For now, the agreement between our different numerical approaches supports the general 
arguments made by Hernquist & Springel (2003), that the SFR density should peak at 2; > 5, 
mostly independent of the details of star formation and feedback. As we have demonstrated 
here, such an early peak in the cosmic star formation history yields that clearly exceeds 
current observational estimates and the results of semi-analytic models, suggesting that most 
of the stars in the universe at 2; > 3 are "missing." 

We thank Rachel Somerville, Carlton Baugh, Shaun Cole, and Nicola Menci for provid- 
ing the result of their semi-analytic model. This work was supported in part by NSF grants 
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for Astrophysics. The TVD simulations were performed at the National Center for Super- 
computing Applications (NCSA). 
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Tabic 1. Simulations 



Run Lbox 


[h-^ Mpc] 


-^mesh/ptcl 


moM [h ^Mq] mgas [h ^M©] 


M [h-^ kpc] 


TVD: N864L11" 


11.0 


864^ 


1.1 X 10^ 2.7 X 10^ 


12.7 


TVD: N864L22" 


22.0 


864=^ 


8.9 X 10^ 2.2 X 10^ 


25.5 


TVD: N768L25'' 


25.0 


768^ 


2.0 X 10^ 3.4 X 10^ 


32.6 


SPH: Q5^ 


10.0 


3243 


2.1 X 10^ 3.3 X 10^ 


1.2 


SPH: 05^= 


33.75 


324^ 


8.2 X 10^ 1.3 X 10^ 


4.2 


SPH: G6^ 


100.0 


486^ 


6.3 X 10« 9.7 X 10^ 


5.3 



Note. — Parameters of some of the simulations on which this study is based. The quan- 
tities hsted are as follows: Lbox is the simulation box-size, A^mesh/ptci is the number of the 
hydrodynamic mesh points for TVD or the number of gas particles for SPH, moM is the dark 
matter particle mass, mgas is the mass of the baryonic fluid elements in a grid cell for TVD or 

the masses of the gas particles in the SPH simulations. Note that the TVD uses 432^ (384^) 
dark matter particles for N864 (N768) runs. A£ is the size of the resolution element (cell 
size in TVD and gravitational softening length in SPH in comoving coordinates; for proper 
distances, divide by 1 + z). The upper indices on the run names correspond to the following 
sets of cosmological parameters: (I^m, ^a, ^6, ^, n, ag) = (0.29, 0.71, 0.047, 0.7, 1.0, 0.85) for (a), 
(0.3, 0.7, 0.035, 0.67, 1.0, 0.9) for (b), and (0.3, 0.7, 0.04, 0.7, 1.0, 0.9) for (c). 
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Fig. 1. — Development of stellar mass density as a function of redshift. Panel (a): The 
black solid line is from the most recent TVD N864L22 run (consistent with the not shown 
N864L11 run), and the red short-dashed line is from the TVD N768L25 run which was used 
in Nagamine et al. (2001a). The blue long-dashed line is the analytical model of Hernquist & 
Springel (2003), and the cyan dot-long-dashed line is the result from the SPH 'G6'-run. The 
observational data points are, from low to high redshift. Cole et al. (2001, black filled circle 
at z = 0), Brinchmann & Ellis (2000, cyan open squares aX z = 0.5 — 1), Cohen (2002, black 
open triangles at z = 0.5 — 1), Glazebrook et al. (2004, black filled triangles at 2; = 1 — 2), 
Rudnick et al. (2003, magenta filled squares at z = — 3), and Fontana et al. (2003, blue open 
crosses at z = 0.5 — 3). The four green boxes are from Dickinson et al. (2003a) which show 
the range of systematic uncertainty introduced by varying the metallicity and star formation 
histories of the mass-fitting model they used to derive the stellar mass density. All the data 
points are normalized to the local estimate by Cole et al. (2001) following Rudnick et al. 
(2003). Panel (h): Observational data points and the H&S model are the same as panel (a). 
The results of two semi-analytic models by Somerville et al. (2001, black long-short-dashed 
line), Granato et al. (2000, magenta short-dashed line; GALFORM), and Menci et al. (2002, 
red dot-long-dashed line) are shown. The model of Somerville et al. (2001) is the 'accelerated 
quiescent' model which was used in Somerville et al. (2004), and the model of Menci et al. 
(2002) is the 'no-burst' model that was used in Poll et al. (2003) and Fontana et al. (2003) 
for comparison with their data. 
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Fig. 2. — Star formation rate density as a function of redshift. Panel (a): The hnes corre- 
spond to TVD N864L22 (sohd blue hue), TVD N864L11 (red short-dashed hue; stopped at 
z = 3 owing to its small box-size), TVD N768L25 (black dot long-dashed line), SPH 'D5'- 
run (magenta dot long-dashed line), SPH 'G6'-run (cyan short-dash long-dash line), and the 
model of Hernquist & Springel (2003, green long-dashed line) as given by equation (1). The 
simulation results shown in this figure are extracted directly from the runs, without any 
further adjustments. The result of the SPH Q5-run is not shown here because it follows the 
Hernquist & Springel (2003) model line very closely at 2; > 3. The sources of the (dust cor- 
rected) observational data points are: Lilly et al. (1996, filled circles), Ferguson et al. (1996); 
Madau (1997, filled squares), Connolly et al. (1997, open squares), Sawicki et al. (1997, filled 
triangles), Treyer et al. (1998, open cross), Pascarelle et al. (1998, open circles), Cowie et al. 
(1999, open diamond), Steidel et al. (1999, open stars at z = 3,4), Iwata et al. (2003, open 
pentagon at 2; = 5), Ouchi et al. (2003b, open inverted triangle at z = 4, 5), and Giavalisco 
et al. (2004, open triangles at 2; = 3 — 6). Panel (h): H&S model (green long-dashed line), 
semi-analytic models of Somerville et al. (2001, blue solid line), Granato et al. (2000, black 
short-dashed line), and Menci et al. (2002) are compared with observations. The model of 
Somerville et al. (2001) is the 'accelerated quiescent' model which was used in Somerville 
et al. (2004), and the model of Menci et al. (2002) is the 'no-burst' model that was used in 
Poll et al. (2003) and Fontana et al. (2003) for comparison with their data. Both simulations 
and H&S model has a peak of SFR at 2 = 5 — 7, whereas the semi-analytic models have a 
peak at 2; = 2 — 4. 



